Baltic Astronomy, vol. 8, 535-574, 1999. 



GALAXY SURFACE PHOTOMETRY 



Bo Milvang-Jensen^'^ and Inger J0rgensen^ 

^ Copenhagen University Observatory, 2100 Copenhagen 0, Denmark, 

milvangOastro .ku.dk 
^ School of Physics and Astronomy, University of Nottingham, University 

Park, NC7 2RD Nottingham, UK (postal address for BMJ) 
^ Cemini Observatory, 670 N. Aohoku PL, Hilo, Hawaii 96720, USA, 

i j orgensenOgemini . edu 

Received March 3, 2000 



Abstract. We describe galaxy surface photometry based on fitting 
elfipses to the isophotes of the galaxies. Example galaxies with differ- 
ent isophotal shapes are used to illustrate the process, including how 
the deviations from elliptical isophotes are quantified using Fourier 
expansions. We show how the definitions of the Fourier coefficients 
employed by different authors are linked. As examples of applications 
of surface photometry we discuss the determination of the relative 
disk luminosities and the inclinations for E and SO galaxies. We also 
describe the color-magnitude and color-color relations. When using 
both near-infrared and optical photometry, the age-met allicity de- 
generacy may be broken. Finally we discuss the Fundamental Plane 
where surface photometry is combined with spectroscopy. It is shown 
how the FP can be used as a sensitive tool to study galaxy evolution. 

Key words: techniques: photometric - galaxies: photometry - 
galaxies: fundamental parameters - galaxies: elliptical and lentic- 
ular, cD - galaxies: stellar content - galaxies: evolution 

1. INTRODUCTION 

Surface photometry of galaxies is a technique to quantitatively 
describe the light distribution of the galaxies, as recorded in 2- 
dimensional images. This paper focuses on the techniques used to 
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derive the surface photometry and presents some examples of scien- 
tific apphcations. The paper is a summary of the lectures given at 
the summer school, and is not intended as a complete review of the 
topic. Surface photometry of galaxies has previously been reviewed 
by Kormendy & Djorgovski (1989) and Okamura (1988). 

The techniques and the software used for performing surface 
photometry of galaxies are described in Section 2. Surface photom- 
etry has many applications. In Section 3 we give an example of 
such an application, namely the determination of disk luminosities 
and inclinations of E and SO galaxies. This determination is based 
on surface photometry - ellipticities and deviations from elliptical 
shape - alone. 

From surface photometry in several passbands we can derive the 
colors of the galaxies. These colors provide information about the 
ages and the metal content of the stellar populations in the galaxies. 
In order to study this, stellar population models are needed. These 
models are the subject of Section 4. In Section 5, we then present 
examples of the color-magnitude and color-color relations of galaxies. 

Although the topic of this summer school is photometry, we will 
nevertheless show an example of the science that can be carried out 
when surface photometry is combined with spectroscopy. The chosen 
example is the relation known as the Fundamental Plane for E and 
SO galaxies, which we discuss in Section 6. 

We end this description of surface photometry with a few sug- 
gestions for future projects that can be carried out based on surface 
photometry only, see Section 7. 

Throughout the paper we use Hq = 50kms~^ Mpc~^. 

2. SURFACE PHOTOMETRY 

Surface photometry is used to study extended objects, such as 
galaxies, as opposed to point sources, such as stars. From an image 
of a point source, only its total magnitude can be derived. From an 
image of a galaxy, it is possible to determine a number of quantities. 
Some of these quantities are derived from surface photometry, such 
as how the intensity and ellipticity vary with radius. Other quanti- 
ties are determined by other methods - the morphological type, for 
example, is determined from visual inspection of the image. There 
also exist schemes to do automated morphological classification, e.g. 
Abraham et al. (1994, 1996); Naim et al. (1995). 
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In this context, we are talking about galaxies where the individ- 
ual stars cannot be resolved in the images. This is for example the 
case for the Hydral (Abell 1060) cluster, which is a nearby cluster 
at a distance of ~ 80 Mpc. At that distance, even the Hubble Space 
Telescope (HST) cannot resolve individual stars in the galaxies. 

Photometry of crowded stellar fields, such as globular clusters, 
is not considered surface photometry. Surface photometry is used 
when the magnitudes of the individual stars (typically in a galaxy) 
cannot be measured, but only their smooth integrated light. 

2.1 Ellipse fitting 

Surface photometry of galaxies is usually done by fitting ellipses 
to the isophotes. This choice is motivated by that fact that the 
isophotes of galaxies are not far from ellipses. This is especially the 
case for elliptical (E) and lenticular (SO) galaxies. In this paper 
we will concentrate on E and SO galaxies. We will also limit our 
discussion to data obtained with CCDs (Charge- Coupled Devices). 

There exist several software packages for deriving surface pho- 
tometry. For the exercise related to the lectures given at the summer 
school we have chosen the ELLIPSE task in the ISOPHOTE pack- 
age (see Busko 1996). This task is based on the ellipse fitting algo- 
rithm used in the GASP package by Cawson (1983; see also Davis 
et al. 1985), the code for which was later rewritten by Jedrzejewski 
(1987). ISOPHOTE is a part of the external IRAF* package STS- 
DAS**. We chose this software for the exercise since it is publicly 
available, and because it includes some documentation. The exam- 
ples of surface photometry presented in this section are also based 
on the ISOPHOTE package. Our aim is to describe the basic prin- 
ciples of surface photometry based on ellipse fitting, and the choice 
of software package is not critical for that purpose. A comparison of 
the results obtained with ELLIPSE with those obtained with other 
software packages is beyond the scope of this paper. 

* IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. See also [http : / / iraf . noao . edu7 

** STSDAS is distributed by the Space Telescope Science Institute, 
operated by AURA, Inc., under NASA contract NAS 5-26555. See 
also [http : //ra. stsci . edu/STSDAS . html 
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To illustrate how surface photometry on E and SO galaxies is 
derived we have chosen three example galaxies, see Figure 1 and Ta- 
ble 1. We will refer to these galaxies as the 'pure E', the 'SO' and the 
'boxy E', respectively. The name 'pure E' refers to the fact that the 
isophotes of this galaxy are almost perfect ellipses. For the 'boxy E', 
the isophotes are box shaped. The pure E has been morphologically 
classified as E3/S0. The boxy E has been morphologically classified 
as SB(rs)0(0) (barred SO with rings). For our purpose of illustrating 
surface photometry these morphological types are not important. 



Table 1. The three example galaxies 



Description 


Name 


Morph. type 




re 


Pure E 


R347 / IC2597 


E3/S0 


—23.35 mag 


8.9 kpc 


SO 


R338 


S0(5) 


-20.74 mag 


1.9 kpc 


Boxy E 


R245 


SB(rs)0(0) 


-20.94 mag 


4.8 kpc 



Note: Name and morphological type are from Richter (1989). To- 
tal absolute Gunn r magnitude (Mj-^) and effective radius (r^) are from 
Milvang- Jensen & J0rgensen (2000), based on Hq — 50kms~^ Mpc~^. 



The example galaxies are all members of the Hydral (Abell 1060) 
cluster. The images used are 300 sec Gunn r exposures (Aes = 
6550 A, Thuan & Gunn 1976) obtained with the Danish 1.5-meter 
Telescope, LaSilla, Chile. The spatial scales, as well as the grey- 
scales of the images of the three galaxies in the figures are identical, 
allowing a direct visual comparison between the three galaxies. The 
intensity scaling in the images shown on Figure 1 is logarithmic. We 
have given lengths in kpc rather than arcsec. To calculate lengths 
in arcsec, use log(£/arcsec) = log(£/kpc) -1-0.407. The corresponding 
distance modulus for the cluster is (m — M) = 34.60 mag. 

In doing the surface photometry, the first step is to identify 
the other objects in the image (other galaxies, stars, and cosmic ray 
events) and mask (fiag) these. This is shown in Figure 2. The masked 
pixels are not used in the surface photometry of the galaxy. The 
masking shown in Figure 2 uses squares for the masking of objects. 
A masking using circles is of course also possible, and that would 
mask fewer 'uncontaminated' pixels. 

Once the masking is done, we need only provide the surface 
photometry task the approximate location of the center of the galaxy. 
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Fig. 1. The three example galaxies - pure E (top), SO (left), boxy E 
(right). The width of the montage is 32 kpc. North is down and east is to 
the right. 
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The task then fits eUipses to the intensities in the image. This is 
done at a number of discrete radii, as also shown in Figure 2. For 
the ELLIPSE task, these discrete radii are specified by the rule that 
the different semi-major axis lengths a are spaced by a factor of 1.1. 

The center {x,y) and the shape (ellipticity e and position angle) 
of the ellipses are kept as free parameters in the fit for semi-major 
axes at which the signal-to-noise ratio is sufficiently high. For larger 
semi-major axes, where the signal-to-noise ratio is lower, the center 
and shape of the ellipses are fixed. 

We get two types of output from the ellipse fit. One type is the 
residual image, which is the difference between the original image 
and the model image based on the best-fitting ellipses. The residual 
images for the three example galaxies are shown in Figure 3. 

For the pure E, the residuals are fairly small. For the SO in 
particular, and for the boxy E, the residuals are larger. The pure 
E is much brighter than the other two galaxies, so in relative terms 
the residuals for the pure E are much smaller than for the other 
two galaxies. For the pure E, very little structure is seen in the 
residual image. The residual images of the SO and the boxy E show 
clear structures. In Section 2.2 we will discuss how to quantify these 
structures. 

The other type of output from the ellipse fit is the radial profiles 
of a number of quantities. I.e., for each ellipse we get the inten- 
sity, center, ellipticity, position angle, and measures of the deviations 
from perfect elliptical isophotes (see Section 2.2), as well as the un- 
certainties for all these quantities. In Figure 4 we show the radial 
profiles of intensity, ellipticity and position angle for the three ex- 
ample galaxies. The position angles shown in Figure 4 are measured 
counter-clockwise from the y-axis of the images. The ELLIPSE task 
adopts this definition of position angles. The standard (astronom- 
ical) definition of position angles is from north through east. Our 
images have north down and east to the right, and thus need to be 
rotated by 180° to have the y-axis poiting towards north. However, 
position angles are only unique to within 180° since the major axis 
of a galaxy does not have a direction (as opposed to a coordinate 
axis). Therefore, the position angles shown in Figure 4 are expressed 
in the standard way. It also follows, for example, that the position 
angle of the inner isophotes of the boxy galaxy could be said to be 
150° as well as -30°. 

In Figure 4, we have plotted the different quantities against the 
logarithm of the equivalent radius r. The equivalent radius is defined 
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Fig. 3. The three example galaxies: residual images (i.e. the original 
image minus the model image). The intensity scaling is linear, and the 
cuts are symmetrical around zero. Black represents positive residuals. 
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Fig. 4. Radial profiles of intensity, ellipticity and position angle for 
the three example galaxies. The plotted range in equivalent radii r is 
0.3-40 kpc (0.8"-102")- 



by r = Vab, where a and b arc the semi-major and semi-minor axis of 
the ellipse, respectively. A circle with radius r has an area equivalent 
to an ellipse with axes a and b, hence the name equivalent radius. 
The program used for the illustrations in this section, ELLIPSE, 
uses a to characterize the size of the ellipses. However, since the 
ellipticity e of the ellipse is defined as e = 1 — 6/a, it follows that r 
can be calculated from a and e as r = ay/1 — e. 

From Figure 4 it is seen how the ellipticity and position angle 
are free parameters until a certain radius were their values are fixed. 
The ellipticity is seen to vary with radius for all three galaxies. The 
position angle for the boxy galaxy is seen to vary rapidly in the outer 
parts, which is also seen in Figure 2. This behavior is known as an 
isophote twist. 

2.2 Quantifying the deviations from elliptical shapes 

As we saw from the residual images in Figure 3, the isophotes 
of E and SO galaxies are not always perfectly elliptical. We wish to 
quantify these deviations from elliptical shapes. To illustrate how 
this is done, we have chosen an example ellipse for each of the three 
example galaxies. The three example ellipses are shown in Figure 5, 
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6 and 7 - they are overlayed both on the original and on the residual 
images of the example galaxies. 

In Figures, 6 and 7, we also plot the quantity 



AI = (1) 

norm — I dl \ ' v / 



dJ 
dr I 




Data 



Data(3): intensity not normalized 



100 200 300 

6 — Angle CCW from major axis [deg] 
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2800 
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2400 
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d — Angle CCW from major axis [deg] 



Fig. 5. Top: Original and residual images of the pure E galaxy, with 
a single ellipse at a = 2.6 kpc (r = 2.3 kpc) shown. The images are 13 
kpc on the side. The 'Data' plot shows A/norm versus 9 (see text) along 
the shown ellipse. The 'Data(2)' plot shows the same, except that the 
intensity / is used rather than A/norm • Little structure is seen, consistent 
with the fact that all the Fourier coefficients are close to zero, see Table 2. 
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Fig. 6. Top: Original and residual images of the SO galaxy, with a 
single ellipse at a = 2.2kpc (r = 1.4kpc) shown. The images are 13 kpc 
on the side. The 'Data' plot shows A/norm versus 9 along the shown 
ellipse. Substantial structure is seen. The starting point of the curve 
{6 = 0°) corresponds to the apogee of the ellipse that is closest to the 
top of the figure. At that point the residual is positive (black). The angle 
9 is measured counter-clockwise. The two most dominant Fourier modes, 
the 4th and 6th order cosine terms, cf. Table 2, are illustrated in the three 
other plots. This example galaxy is a disky galaxy and as such has a 
positive C4 (4th order cosine coefficient). 



546 B. Milvang- Jensen & I. J0rgensen 




Data 4lh order cosine term 




B [deg] S [deg] 

Fig. 7. Top: Original and residual images of the boxy galaxy, with 
a single ellipse at a = 3.8 kpc (r = 2.8 kpc) shown. The images are 13 

kpc on the side. The 'Data' plot shows A/norm versus 9 along the shown 
ellipse. Substantial structure is seen. The starting point of the curve 
= 0°) corresponds to the apogee of the ellipse that is closest to the top 
of the figure. At that point the residual is negative (white). The angle 
9 is measured counter-clockwise. The two most dominant Fourier modes, 
the 4th order cosine and sine terms, cf. Table 2, are illustrated in the three 
other plots. Note how a boxy galaxy has a negative C4 (4th order cosine 
coefficient) . 
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As can be seen, A/norm is the deviation in intensity / from the 
mean intensity at the given eUipse, Iq, normahzed by the equivalent 
radius r, and by the absolute value of the intensity gradient, |d//dr|. 
For comparison, we also show the intensity / in Figure 5. The reason 
for choosing this particular normalization will be described in the 
following. One advantage of using the quantity A/norm in the plots 
is that allows a direct comparison between the plots for the three 
example galaxies. 

To quantify how the intensity deviates from being constant along 
the fitted ellipse, the following Fourier series is fitted to the intensity 

m 

N 

/(^) = Jo + J] [A^ sin(n^) + 5^ cos(n^)] . (2) 

n=l 

N is the highest order fitted. 9 is the angle measured counter- 
clockwise from the major axis of the ellipse. The diff'erent Fourier 
modes, e.g. cos 4^, will be discussed below. 

We are more interested in how the isophote deviates (in the 
radial direction) from the fitted ellipse. Let Riso{0) denote the dis- 
tance from the center of the ellipse to the isophote, and let i?eii(^) 
denote the distance from the center of the ellipse to the ellipse itself. 
For a perfectly elliptical isophote, the difference between Riso{0) and 
Re\i{0) would be zero for all values of 6. We can Fourier expand the 
difference as 

N 

AR{e) = Ri,o{d)-Reii{0) = J2[Ksm{ne) + B'^cos{ne)] . (3) 

n=l 

AR{9) is the radial deviation of the isophote from elliptical shape. 
The relative deviation is more interesting, so we take AR{9) relative 
to the size of the ellipse, given by the equivalent radius r. This rela- 
tive radial deviation of the isophote from elliptical shape, AR{9)/r, 
is described by the Fourier coefficients A'^/r and B'^/r. We will al- 
locate new symbols for these quantities: s„ = A'^/r and c„ = B'^/r. 

The ELLIPSE task calculates An and but not A'^ and B'^,. 
However, we are able to link the two sets of coefficients. Consider a 
Taylor expansion to first order of /(i?) around Rq 



I{R)^I{Ro) + -^{R-Ro) 



(4) 
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Let -R be a point on the isophote and let Rq be a point on the 
eUipse. I{R) is constant since R is on the isophote. The intensity 
on the elhpse I{Ro) is not necessarily constant. The intensity on 
the ellipse is also given by I (9) (Equation 2). We can identify the 
difference {R — Rq) with AR{9) (Equation 3). A suitable mean value 
of the gradient dl/dR is d//dr. (The 'effective intensity gradient' 
d//dr can be calculated as the difference in intensity divided by the 
difference in equivalent radius for two adjacent ellipses.) By inserting 
Equation (2) and (3) in Equation (4) we find the following relations 
hold for all n 

I dr I I dr I 

where we have used — d//dr = |d//dr| since the intensity gradient is 
negative. From the definitions of s^i and c„ given above, we finally 
get 

s - ^" c - (6) 

' ■ I dr I ' ■ I dr I 

These definitions of the Fourier coefficients Sn and c„ are used 
in the literature by e.g. Franx, lUingworth & Heckman (1989b); 
J0rgensen, Franx & Kjaergaard (1992); J0rgensen & Franx (1994); 
and J0rgensen, Franx & Kjaergaard (1995). 

Slightly different definitions are also in use. Some authors, 
e.g. Bender & MoUenhoff (1987); Bender, Dobereiner & Mollenhoff 
(1988); Bender et al. (1989); and Nieto & Bender (1989), use 04/0 
for the 4th order cosine coefficient (still for the radial deviation). In 
our notation, this is equal to B'^/a. The only difference between C4 
and a^/a is that C4 = B^r is taken relative to the equivalent radius 
r, whereas 04/0 = -84/0 is taken relative to the semi-major axis a. 
Thus, the two are related by 



a^/a — C4^r/a — C4VI — e = Ci^\/hja , (7) 

where we have used the known relations between r, a, h and e. For 
an apparently round galaxy (i.e. for e = 0), a^ja is equal to C4. 

Yet another definition is used by e.g. Peletier et al. (1990). These 
authors expand the intensity along the ellipse as 

AT 

I{Q) = Jo (1 + X! \.^n sin(n^) + Cn cos(n^)]) (8) 

n=l 
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(note the upper case C„). When comparing with Equation (2) it is 
seen that loSn = An and loCn = B^- This means, for example, that 
C4 and C4 are related as 



C4 = C4 



r 


dl 




dlog/ 




dr 


= C4 • 


dlogr 



(9) 



It also follows that a^/a is related to C4 by 



a^/a = C4 • 



d log/ 



dlogr 



(10) 



a relation used by Faber et al. (1997) (but note the upper case C4). 



Table 2. Fourier coefficients for the three example galaxies at 
the particular ellipses shown in Figure 5, 6 and 7. 
The semi-major axis a for these ellipses is given. 



Coeff. 




a = 


pure E 
= 2.6 kpc 




so 

a = 2.2 kpc 




a = 


boxy E 
= 3.8 kpc 


C3 





001 


± 0.002 


-0 


003 ±0.013 





002 


± 0.008 


S3 


-0 


001 


± 0.002 


-0 


000 ± 0.014 


-0 


002 


± 0.008 


C4 





000 


± 0.002 





065 ±0.013 


-0 


058 


± 0.006 


S4 


-0 


004 


zb 0.002 





001 ± 0.008 


-0 


032 


±0.005 


C5 





000 


± 0.002 


-0 


003 ± 0.008 





002 


± 0.004 


S5 





001 


± 0.002 





007 ± 0.008 


-0 


006 


± 0.004 


C6 





003 


±0.001 





037 ±0.006 


-0 


018 


± 0.004 


S6 


-0 


002 


±0.001 


-0 


000 ± 0.006 





012 


± 0.004 


C7 


-0 


000 


±0.001 


-0 


001 ± 0.006 





010 


±0.003 


S7 


-0 


002 


±0.001 


-0 


000 ± 0.006 


-0 


005 


± 0.003 


C8 





004 


±0.001 





021 ± 0.005 


-0 


004 


± 0.003 


S8 


-0 


001 


±0.001 





003 ± 0.005 





005 


± 0.003 



The fitted values of the Fourier coefficients Cn and Sn for the 
example ellipses for the three example galaxies are listed in Table 2. 
For the SO galaxy, it is seen that the two numerically largest coeffi- 
cients are C4 = 0.065 and cq — 0.037. For the boxy galaxy the two 
numerically largest coefficients are C4 = —0.058 and S4 = —0.032. 

From the definition of A/norm (Equation 1), and from Equation 
(3) and (4) it is seen that A/norm = Ai?(6')/r, i.e. A/norm measures 
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Fig. 8. Radial profiles of the Fourier coefficients C4, S4, Cg and Sq 
for the three example galaxies. The plotted range in equivalent radii r is 
0.3-40 kpc (0.8"-102")- Points with uncertainty larger than 0.025 are not 
plotted. 



the relative radial deviation of the isophote from elliptical shape. 
This is also the quantity measured by Sn and Cn (indeed, they are 
the Fourier coefficients of AR{6)/r), so plots of the Fourier modes 
Sn sin(n6') and Cn cos{n9) can be directly compared with the A/norm 
versus 9 plot. In Figure 6 and 7 we show the two most important 
Fourier modes and their sum for the SO galaxy and the boxy E galaxy, 
respectively. It is seen in both cases how the two most important 
modes account for most of the structure. 

An ellipse can be described by the first and second order Fourier 
coefficients. Since the Fourier expansion is along the best-fitting 
ellipses, the first and second order Fourier coefficients obtained there 
will be zero. 
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The cosine modes are symmetrical around the major axis {9 = 0° 
or 180°), whereas the sine modes are not, see Figure?. The odd co- 
sine modes (e.g. 3rd order) are not symmetrical around the minor 
axis (see Figure la in Peletier et al. 1990). The mode that is dom- 
inating for most E and SO galaxies is the 4th order cosine. Note 
from Figure 6 and 7 how C4 is an indicator of disky (04 > 0) or boxy 
(c4 < 0) isophotes aligned with the major axis (e.g. Carter 1987; 
Bender et al. 1989; Peletier et al. 1990). 

The radial profiles of the 4th and 6th order Fourier coefficients 
are shown in Figure 8. 

2.2.1 The output from the ELLIPSE task 

The radial profiles determined by the ELLIPSE task are output 
to an STSDAS table. The most important columns are listed in 
Table 3, along with their meaning in our notation. 

The equivalent radius is not in the table, but it can be calculated 
as SMA*sqrt (1-ELLIP). The 5th to 8th order Fourier coefficients s„ 
and Cn can be calculated as e.g. cq = BI6/ (SMA*abs (GRAD) ) . This 
follows from Equation (6) since ad//da = rd//dr. Note that the 
GRAD column is d//da (I. Busko 2000, private communication). 

Table 3. Output columns from ELLIPSE 



Column name Content in our notation 



SMA 


a 


INTENS 


I 


ELLIP 


e 


PA 


PA (position angle) 


XO 


X (center of ellipse) 


YO 


y (center of ellipse) 


GRAD 


d//da (not d//dr) 


An 


Sn, n = 3, 4 


Bn 


Cn, n = 3, 4 


Aln 


A^, n = 5, 6, 7, 8^ 


Bin 


Bn, n = 5, 6, 7, 8^ 



^ Only calculated with the following option set: 
harmonics="5 6 7 8". 
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2. 3 Determination of magnitudes 

The intensity / (in ADU) contains the signal from the galaxy 
plus the sky background. The sky background level can be deter- 
mined in several ways. One way is to identity empty regions in the 
image and measure the level there. If the galaxy fills most of the im- 
age, this can be difficult. Another way is to fit a suitable analytical 
expression to outer part of galaxy plus sky intensity profile obtained 
from the surface photometry. This method has been used by e.g. 
J0rgenscn et al. (1992), fitting Igaiaxy-i-sky (r) = Igky + -^gaiaxy.o • r'"', 
with a = 2 or 3. 

Magnitudes can be calculated from the sky subtracted intensity. 
This can either be integrated magnitudes within a certain aperture 
(elliptical or circular), or the surface brightness at a given ellipse, 
n{r). By knowing the pixel scale of the CCD (in arcsec/pixel), //(r) 
can be expressed in units of mag/arcsec^. With the use of observed 
standard stars, the magnitudes and surface brightnesses can be trans- 
formed to a standard photometric system. 

2.4 Global parameters 

The surface photometry has produced radial profiles of a number 
of quantities. It is desirable to condense these radial profiles to a few 
characteristic numbers, the global parameters. 

2.4.1 Effective parameters 

Elliptical galaxies have surface brightness profiles that arc well 
approximated by the r^/^ law (de Vaucouleurs 1948). By fitting 
the aperture magnitudes to an r^/*^ growth curve, the following two 
parameters can be derived: 

• re'- Effective radius, in arcsec 

• (/u)e: Mean surface brightness within re, in mag/arcsec^ 
The seeing needs to be taken into accoimt (Saglia et al. 1993). 

For galaxies with perfect r^/^ profiles, the effective radius r^ is 
the half-light radius, i.e. the radius that encloses half of the light from 
the galaxy. Spiral galaxies are better described by an exponential 
surface brightness profile than by an r^/^ profile. 

We can express the mean surface brightness in units of Lq/pc^, 
where Lq is the luminosity of the Sun in the given passband 
(e.g. Gunn r). We will call this quantity (/)e. The relation is 
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log(/)e = — 0.4((/x)e — k), where the constant k is given by /c = 
Mq + 51og(206265pc/10pc). As is seen, the calculation does not 
involve the distance to the galaxy, but only the absolute magnitude 
of the Sun in the given passband. The calculation of in kpc from 
Te in arcsec, however, does involve the distance to the galaxy. 

With Te as the half-light radius, it follows that the total lumi- 
nosity is given by L = 27r {I)^r^. 

2.4.2 Global Fourier parameters 

As is seen from Figure 8, also the Fourier parameters vary with 
radius. One way of getting a characteristic value of e.g. the C4(r) 
profile is to take the extremum value. In case the profile does not 
have a clear extremum, we can take the value at the effective radius. 
We will use the symbol C4 for this characteristic value of C4(r). 

Another way to get a global Fourier coefficient is to calculate an 
intensity weighted mean value as 



where rmin is the radius where seeing effects are no longer important, 
and Tmax is the radius where the Fourier coefficients can no longer 
reliably be determined (see J0rgensen & Franx 1994). 

2.4.3 Global ellipticities and colors 

As global ellipticity can be taken the extremum of the e(r) pro- 
file, the value of e(r) at the effective radius, or the value of e(r) at a 
certain isophote level, e.g. /j, = 21.85'^/arcsec^ in Gunn r as used by 
J0rgensen & Franx (1994). 

As global color, the color within the effective radius can be used. 
By color is meant the difference between the magnitudes in two dif- 
ferent passbands, such as (B— V). The color is always calculated as 
the magnitude in the passband with the shortest effective wavelength 
minus the magnitude in the passband with the longest effective wave- 
length. Thus, a large value of the color means that the galaxy is red, 
and a small value means that it is blue. Another global parameter 
related to the color is the color gradient, defined as Acolor/A logr, 
i.e. as the slope of the color versus logr plot. 




(11) 
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3. QUANTITATIVE MORPHOLOGY FOR E AND SO GALAXIES 

CCD surface photometry as described in the previous section 
offers the possibility of carrying out quantitative morphologic studies 
of E and SO galaxies. The traditional classification of these galaxies 
distinguishes between E and SO galaxies based on the prcsense of a 
disk (SO galaxies) or no disk (E galaxies). E and SO galaxies are 
considered to belong to two separate classes of galaxies. 

In the following we summarize the methods and results presented 
by J0rgensen & Franx (1994, hereafter JF94). The results lead to 
the conclusion that E and SO galaxies fainter than an absolute blue 
magnitude of Mb^ = —22 mag form one class of galaxies with a 
broad and continuous distribution of the relative luminosity of the 
disks. 

3.1. Sample properties and the data 

The sample used in the study by JF94 is a magnitude limited 
sample of galaxies within the central square degree of the Coma 
cluster. The sample is complete to an apparent magnitude in Gunn r 
of 15.3 mag. 171 galaxies are included in the sample. Because the 
sample is well-defined and complete it is possible based on these data 
to draw conclusions regarding the E and SO galaxies as a class. 

CCD photometry was obtained of the full sample. Surface pho- 
tometry for the galaxies was derived using GALPHOT (J0rgensen et 
al. 1992; Franx et al. 1989b). From the surface photometry, global 
parameters were derived. The important parameters used in this 
study are summarized in Table 4. (04) and (ce) are intensity weighted 
mean values of C4 and ce, respectively (see JF94 and Section 2.4.2), 
while C4 represents the extremum the C4 radial profile. The intensity 
weighted parameters are less sensitive to small scale features in the 
radial profiles and are therefore to be preferred for studies of global 
properties. 

3.2. Morphology of the E and SO galaxies 

Figure 9 shows the morphological parameters 621.85, (^4) and (cg) 
as functions of the total magnitudes mx in Gunn r. The spiral galax- 
ies are also shown on these figures, but we will omit them from the 
following discussion. Prom Figure 9 it is clear that the galaxies fainter 
than about 12.7 mag span the full range in morphological parame- 
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Table 4. Surface photometry parameters 



Parameter Description 



rriT 


Total magnitude in Gunn r 


^21. 85 


Ellipticity at fj, — 21.85™/arcsec^ in Gunn r 


C4 


Extremum of the C4(r) radial profile 


(C4) 


Intensity weighted mean value of C4(r) 


ice) 


Intensity weighted mean value of CQ{r) 



ters. The E and SO galaxies cannot be separated into two classes of 
galaxies based on one of these morphological parameters. The dif- 
ference in properties of E and SO galaxies fainter and brighter than 
TTiT = 12.7 mag is striking. This demarcation magnitude corresponds 
to a blue absolute magnitude of Mb^ = —22 mag. 

In Figure 10 we show the distribution of the ellipticities, €21.85, 
of the galaxies. Figure 10b shows the cumulative frequency of 621.85 
for the E galaxies and the SO galaxies separately, as well as the 
cumulative frequency of €21.85 for the E and SO galaxies as one class. 

If the E galaxies and the SO galaxies form two separate classes, 
then we expect that their 621.85 distributions can be modeled inde- 
pendently. The simplest assumption is that the galaxies are ran- 
domly oriented in space and have some simple distribution of the 
intrinsic ellipticities. JF94 attempted to fit the 621.85 distributions 
with intrinsically uniform ellipticity distributions and with intrinsic 
ellipticity distribution that were Gaussians. The fit was done by 
maximizing the probability that the data was drawn from the model 
as reflected by a Kolmogorov-Smirnov (K-S) test (e.g. Press et al. 
1992). The K-S test gives the probability that a data distribution 
is drawn from a model distribution (or another data distribution) 
based on the maximum difference between the cumulative frequen- 
cies of the two distributions, as it is illustrated in Figure 10b. JF94 
found that the 621.85 distribution of the E galaxies could be fitted 
satisfactory with either a uniform or a Gaussian intrinsic distribu- 
tion of the ellipticities. Both of these models resulted in probabilities 
of 84 per cent or larger. However, the 621.85 distribution of the SO 
galaxies could not be fitted satisfactory. Both intrinsic distributions 
have probabilities of only 11 per cent. The 621.85 distribution of the 
SO galaxies lack apparently round galaxies, see Figure 10b. JF94 also 
find that the 621.85 distribution of the E and SO galaxies treated as 
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Fig. 9. Morphological parameters versus the total magnitudes. Open 
boxes - E galaxies with dominating regular C4-profiles; filled boxes - SO 
galaxies with dominating regular C4-profiles; crosses - E galaxies with ir- 
regular or non-dominating C4-profiles; stars - SO galaxies with irregular 
or non-dominating C4-profiles; triangles - spirals. The C4-profiles are con- 
sidered non-dominating if |(c4)| is more than one sigma smaller than the 
absolute value of one of the other intensity weighted mean coefficients. 
Typical measurement errors are given on the panels. Galaxies with un- 
certainty on C4 respectively (04) larger than 0.02 are not plotted. The six 
brightest galaxies have small Fourier coefficients and €21.85 < 0-4. Other 
dependence on ttt-t is not seen. (From JF94.) 



one class can be fitted satisfactory by either a uniform or a Gaussian 
intrinsic distribution of the ellipticities. 

These results show that some (maybe all) of the E galaxies 
fainter than rriT = 12.7 mag must be SO galaxies seen face-on. When 
the galaxies are seen face-on a disk is more difficult to detect (as also 
noted by Rix & White 1990), and the galaxies are mostly classified 
as E galaxies even when a disk is present. 
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Fig. 10. Distributions of the apparent ellipticities. (a) Solid line - E 
galaxies, dashed line - SO galaxies, Dashed-dotted line - spirals, (b) Solid 
lines - E, SO galaxies, and E and SO galaxies together. The six brightest 
E galaxies are excluded. Dotted lines - best fitting uniform distributions. 
Dashed lines - best fitting Gaussian distributions. (From JF94.) 

3.3 The relative disk luminosities 

The next task is to determine the relative disk luminosities, 
L-o/Ltot, of the galaxies. Ld/ L^ot is the luminosity of the disk rela- 
tive to the total luminosity of the galaxy. JF94 showed that two of 
the morphological parameters can be used together to derive Ld / -C'tot 
if a simple model for the disk and the bulge is assumed. 

Figure 11 illustrates this technique for the Coma cluster sample. 
JF94 constructed models consisting of a bulge with an r^/^ profile 
and a disk with an exponential profile. The bulge is assumed to have 
an intrinsic ellipticity of 0.3, while the intrinsic ellipticity of the disk 
is assumed to be 0.85. JF94 tested models for both equal major 
axis of the two components, aeB = «eD) and for OeB = 0.5aeD- The 
results are not significantly different, so here we will concentrate on 
the OeB = «eD Hiodels. JF94 convolved the model images of bulge 
plus disk with a representative seeing and then analyzed the model 
images in the same way as done for the data. Models with relative 
disk luminosities Lj^/Ltot between zero (no disk) and one (all disk) 
were constructed. Further, the inclination was varied between face- 
on (small inclination, cos i = 1) and edge-on (large inclination, cosz 
= 0) in steps of 0.1 in cos i. 

The results in terms of morphological parameters are shown on 
Figure 11. The models reproduce the general variation of C4, (04) and 
(ce) with ellipticity. They also span reasonably well the section of the 
(c4)-(c6) diagram covered by the data. Models of the kind used by 
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Fig. 11. Fourth and 6th order Fourier coefficients plotted against 
elhpticity and versus each other. Data symbols as in Figure 9. Galaxies 
with uncertainty on C4 respectively (04) larger than 0.02 are not plotted. 
Typical measurement errors are given on the panels. The models with 
tteB — OeD are overplotted. The curves are labeled with L-q/L^q^. The 
dashed lines on (b) mark the inclinations where cosi = 0.1, 0.2, 0.3, 0.4, 
and 0.5 (i.e. i = 84°, 78°, 73°, 66°, and 60°), from top to bottom. No 
dashed line is shown for the i — 90° (i.e. edge-on) models, but this line 
would be above the i = 84° dashed line, connecting the end points of the 
solid lines. It is seen that for a given Lj^ /-^tot the highest inclination gives 
the largest (C4). The non-zero coefficients for the pure disk-model are due 
to the inclusion of seeing effects in the models. (From JF94.) 



JF94 cannot reproduce boxy isophotes of the galaxies. However, the 
Coma cluster sample contains only two galaxies fainter than rriT = 
12™7 that have (04) significantly smaller than zero. For galaxies with 
ellipticities larger than 0.3 and (04) larger than 0.007, the models 
are well separated in (04) versus €21.85, see Figure lib. Thus, these 
two parameters can be used to derive the relative disk luminosities, 
Ln/L^oti and the inclinations, z, of the galaxies. 
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Fig. 12. Relative (a) and cumulative (b) frequency for the relative 
disk luminosities. Solid line - determinations from the e2i.85-(c4) dia- 
gram. The distribution has been normalized with the total number of 

galaxies fainter than = 12^7. Dashed line - model prediction for a 
uniform intrinsic distribution, corrected for the incompleteness due to the 
limits enforced on €21.85 and (C4). Dotted lines - model predictions for a 
uniform intrinsic distribution plus a fraction of diskless galaxies. Models 
for fractions of 0.1, 0.2, and 0.5 are shown. A higher fraction of diskless 
galaxies moves the curve for the normalized distribution downwards in 
both (a) and (b). (Prom JF94.) 



Figure 12 shows the resulting distribution of L^/Ltot- It was 
possible to derive LD/-^tot for 52 of the E and SO galaxies in the 
sample. The cumulative frequency shown on Figure 12b is normal- 
ized to the 140 E and SO galaxies fainter than mx — 12.7 mag. 
Overplotted on Figure 12 are models of the distribution of Ld/Lioi. 
These models are uniform distributions with some fraction of diskless 
galaxies added. The best fitting model is a uniform distribution of 
Lo/I'tot between zero and one with an additional 10 per cent diskless 
galaxies. JF94 find that the resulting distributions of inclinations, 
€21.85 and (04) for this model also fit the observed distributions of 
these parameters. 



3.4 Conclusions from JF94 

The results presented by JF94 illustrate the strengths of studies 
of quantitative morphology based on surface photometry for statisti- 
cally well-defined and complete samples. JF94 were able to show that 
the E and SO galaxies in the Coma cluster (fainter than rriT = 12.7 
mag in Gunn r) form one class of galaxies with a broad (most likely 
uniform) distribution of Lj^/Ltot- This result contradicts the tra- 



560 



B. Milvang- Jensen & I. J0rgensen 



ditional classification of these galaxies into two separate classes. It 
also provides constraints for models for morphological evolution of 
galaxies, since the end-result for the Coma cluster represents one 
scenario that the models need to reproduce. 

4. STELLAR POPULATION MODELS 

Stellar population models are tools for interpreting the inte- 
grated light, such as the colors, observed from galaxies. Ideally, we 
want to determine the mix of stars that give rise to the observations. 
This problem is usually under constrained, so it is necessary to make 
some assumptions regarding how the numbers of different types of 
stars are related. Here we will consider so-called single-age singlc- 
metallicity models, also known as single stellar population (SSP) 
models. In these models, all the stars are formed at the same time, 
with distribution in mass given by the chosen initial mass function 
(IMF), and with identical chemical composition. 

SSP models are based on the following ingredients. First, theo- 
retical stellar isochrones are needed. Isochrones are loci in the theo- 
retical HR-diagram (logTefi, logL) for a stellar population of a given 
age and chemical composition. Second, a conversion is needed be- 
tween the theoretical parameters of Tes, L, log (7 (surface gravity) 
and the metallicity [M/H] to the observable parameters such as col- 
ors. This conversion can be either empirical or theoretical. The 
empirical conversion is based on observations of individual stars in 
our Galaxy for which the 'theoretical parameters' can be inferred, 
and the observable parameters measured. The theoretical conver- 
sion is based on model atmospheres and synthetic spectra. Third, 
the IMF has to be specified. 

An example of SSP models are those presented by Vazdekis et 
al. (1996). These models use the isochrones from the Padova group 
(Bertelli et al. 1994). The conversion from the theoretical to the 
observable parameters is empirical. Models are presented for several 
different IMFs. One IMF is a constant below 0.2 Mq, a Salpeter 
(1955) IMF above 0.6 M©, and a spline in the interval 0.2-0.6 M©. 
The models that we use in Section 5 are based on this IMF. 

When the IMF has been specified, the models have only two 
parameters: age and metallicity. The metallicity can be expressed 
either as the mass fraction in heavier elements, Z, or as the metal 
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abundance [M/H] = log{Z/ZQ) (with Zq = 0.02). The models have 
solar abundance ratios, while this may not be the case for all galaxies. 
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Fig. 13. Example of the predictions from tiie Vazdekis et al. (1996) 
models: (B— R^) color as function of age (the x-axis) and metallicity (the 
different symbols). Triangles - [M/H] =0.4; crosses - [M/H] = 0.0; boxes 
- [M/H] = —0.4. It is seen that the stellar populations get redder with 
higher age and higher metallicity. 



As an example, the model predictions for the (B— Rc) color is 
given in Figure 13. It is seen that the color depends both on age and 
metallicity. As will be illustrated in the next section, optical colors 
depend roughly in the same way on age and metallicity. Therefore, in 
a plot of one optical color versus another optical color the model lines 
of constant age will be almost on top of the model lines of constant 
metallicity. The effects of age and metallicity cannot be disentangled 
in such a diagram. This is known as the age-metallicity degeneracy 
(e.g. Worthey 1994). 

It should be noted that real galaxies are not necessarily single 
stellar populations. For example, a galaxy could have experienced a 
second star formation event. Therefore, when SSP model predictions 
are compared with data for real galaxies to determine the age and 
the metallicity, the resulting ages and metallicities are luminosity 
weighted mean values. Further, dust can also cause a stellar pop- 
ulation to appear red. This can be an additional complication in 
determining ages and metallicities. 
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5. COLOR RELATIONS 

The E and SO galaxies follow very well-defined relations between 
tfie optical colors and the total magnitudes. Figure 14 shows the 
optical color-magnitude relation for the Coma cluster. The figure 
includes all objects detected in a field covering the central 75 arcmin 
X 80 arcmin of the cluster. The data were obtained with the Mc- 
Donald Observatory 0.8- meter Telescope equipped with the Prime 
Focus Camera. 

The color-magnitude relation is well-defined and has a very low 
scatter for (E and SO) galaxies brighter than about Rc = 17 mag, 
see Figure 14. For galaxies fainter than Rc — 17.5 mag only sparse 
redshift information is available, and many of these faint galaxies 
may be background galaxies. 

The color-magnitude relation is thought to be primarily a result 
of differences in metallicity as a function of luminosity. However, re- 
cent results based on spectroscopy show that both the mean ages and 
the mean metallicities varies for E and SO galaxies at low redshifts 
(e.g. Worthey, Trager & Faber 1995; J0rgensen 1999). Thus, there 
is a need for interpreting the color-magnitude relation within these 
recent results in order to achieve an self-consistent interpretation of 
the spectroscopic and the photometric results. The work by Kauff- 
mann & Chariot (1998) represents one of the only attempts to model 
the color-magnitude relation and relations involving spectroscopic in- 
formation in a consistent manner. In the models by Kauffmann & 
Chariot both age and metallicity varies with the luminosity. 

Color-color relations in the optical are in Figure 15 shown for 
the same sample of objects as shown in Figure 14. The confirmed 
members of the cluster form a tight relation in these two color-color 
relations. This is in agreement with predictions from stellar popu- 
lation models (see Section 4), which predicts the optical colors to 
be degenerate in age and metallicity. This is seen on Figure 15, left 
panel, where the same models as shown on Figure 13 are overplotted. 
The lines of constant age fall right on top of the lines of constant 
metallicity, forming a single line along the ridge of the location of 
the E and SO galaxies. This means that the optical colors alone can- 
not be used to derive the mean ages and the mean metallicities. A 
younger age will lead to bluer colors, but a similar change in colors 
can be caused by a lower metallicity. 

While the optical color-color diagrams arc degenerate in age and 
metallicity, the combination of one optical color and one optical- 
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Fig. 14. Color-magnitude relation for the central 75 arcmin X 80 
arcmin of the Coma cluster. Both stars and galaxies are included on 
this figure, a total of 15370 objects. The lack of objects in the upper 
right hand corner is due to the combination of the magnitude limits in 
B and -Rc- Small black points - stars; black [blue] boxes - spirals and 
irregulars, confirmed members; dark grey [red] boxes - E and SO galax- 
ies, confirmed members; light grey [green] boxes - unclassified confirmed 
members; small grey [orange] points - other galaxies in the field, members 
and non-members. The photometric data are from J0rgensen (2000). The 
redshift data are from J0rgensen & Hill (2000). 



infrared color may be used to break the degeneracy. An example of 
this is shown in Figure 16 for a smaU sample of E and SO galaxies in 
the Coma cluster. The optical color (U— B) and the optical-infrared 
color ( V— K) are both sensitive to both age differences and metalhcity 
differences. However, (U— B) is more sensitive to the age differences 
than to metallicity differences, while the opposite is the case for 
(V— K). The stellar population models in the near infrared (near- 
IR), JHK, in this case the K-band, are stiU rather uncertain, but 
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Fig. 15. Optical color-color relations for the Coma cluster. Symbols 
as in Figure 14. Solid line on the left panel - SPP models from Vazdekis 
et al. (1996), see text. 

the technique is promising for studies of faint high redshift galaxies 
for which spectroscopy would come at a very high cost of telescope 
time at 8-meter class telescopes. The apparently rather low mean 
age of the Coma cluster galaxies as estimated from Figure 16 is in 
fact in agreement with results based on spectroscopic information 
(see J0rgensen 1999). 




6. THE FUNDAMENTAL PLANE 

The Fundamental Plane (FP) is a relation that combines surface 
photometry with spectroscopy. We will discuss this relation both at 
low and at high redshift. 

The FP (Djorgovski & Davis 1987; Dressier et al. 1987) is the 
relation 

logre = alogcr + /31og(/)e + 7 , (12) 

where a is the line-of-sight stellar velocity dispersion for the galaxy 
in question. In other words, the measured values of logrg, log(/)e 
and logo" for a sample of E and SO galaxies do not populate this 
3-parameter space evenly, but are limited to a thin plane. 

The velocity dispersion a is determined from spectroscopy, see 
Figure 17. The absorption lines in galaxies are broadened due to the 
internal motions of the stars in the galaxy. To determine how much 
the stars are moving, it is necessary to know what the spectrum of 
the galaxy would be if all the stars were at rest with respect to each 
other. This is approximated by the spectrum of af K giant star. This 
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Fig. 16. Optical-infrared color-color relation for the Coma cluster. 
The data are from J0rgensen (2000) and Mobasher et al. (1999). The 
lines represent stellar population models by Vazdekis et al. (1996). Solid 

lines - metallicities of [Fe/H] = — 0.4, 0.0 and 0.4; the lowest metallicity 
leads to the smallest (V— K). Dashed lines - ages of 1, 2, 5, 8, 12, 15, and 
17 Gyr; the largest ages lead to the largest (U— B). 

template star spectrum is broadened by a Gaussian broadening func- 
tion until it matches the galaxy spectrum. The velocity dispersion 
of the galaxy is then the dispersion a of the broadening function. In 
more precise terms, this determination of a can be done using the 
Fourier fitting method (Franx, lUingworth & Heckman 1989a) or the 
Fourier quotient method (Sargent et al. 1977). 

6. 1 The interpretation of the Fundamental Plane 

The physics behind the FP can be illuminated by some simple 
arguments (Djorgovski, de Carvalho & Han 1988; Faber et al. 1987). 
Consider the virial theorem for the stars in the galaxy 



We relate the observable quantities Tg, (J and (/)e to the 'physical' 
quantities (i?), {V^) and luminosity L through 



GM 



(R) 



(13) 
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Fig. 17. Illustration of the effect of the velocity dispersion a. The top 
panel shows a K giant star in our own Galaxy. This star is representative 
of the stellar populations in E or SO galaxies. The two lower panels show 
E or SO galaxies in the Hydral cluster. The Mg b absorption line triplet at 
5177 A (rest frame) is broadened in the galaxy spectra. The instrumental 
resolution is 79 km/s. 

re = kn{R) , a^ = ky{V^), L = ki^{I)^rl , (14) 

The parameters /cr, fcy, and /cl reflect the density structure, kine- 
matical structure, and luminosity structure of the given galaxy. If 
these parameters are constant, the galaxies constitute a homologous 
family. Homology means that structure of small and big galaxies is 
the same. 

Combining Equation (13) and (14) gives 
re = A;s(M/L)-V2(/)-i , = (G/cr/cv/cl)-' • (15) 

For homology ko, will be constant. When this relation is compared 
to the observed FP, 

re = constant • ^l-24±0.07^^^-0.82±0.02 (^g) 

(J0rgenscn, Franx & Kjaergaard 1996, in Gunn r), it is seen that the 
coefl&cients of the FP are not 2 and —1 as expected from homology 
and constant mass-to-light ratios. The product ks{M / L)~^ cannot 
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be constant, but has to be a function of a and (/)e. A non-constant 
ks{M/L)~^ can be explained by a systematic deviation from homo- 
logy (ks varies), or a systematic variation of the M/L ratios, or 
both. When homology is assumed, the observed FP coefficients give 
the relation 

M/Lr oc m0-24±0 03 , (17) 

(J0rgensen et al. 1996). The interpretation of the FP is still a matter 
of debate. The M/L oc interpretation seems to be the most 
favored one, although there is some evidence that non-homology may 
play a role too (e.g. Hjorth & Madsen 1995; Pahre, de Carvalho & 
Djorgovski 1998). 

6.2 The evolution of the Fundamental Plane as a function of redshift 

The Fundamental Plane can be used to study the evolution of 
galaxies as a function of redshift. As explained above, the FP may 
be interpreted as a relation between the masses and the M/L ratios 
of the galaxies. Under the assumption that the masses do not change 
with redshift, e.g. no merging takes place, the evolution of the FP 
zero point with redshift can be interpreted as the evolution of the 
M/L ratios. 

Several authors have studied the FP for clusters at redshifts 
higher than 0.1, see Table 5 for clusters and references. Additional 
studies by Pahre, Djorgovski & de Carvalho (1999) and Kelson et al. 
(1999) are soon to be published in refereed journals. 



Table 5. Fundamental Plane studies of cluster with 2; > 0.1 



Cluster 


z 


A^gal 


Reference 


A2218 


0.18 


9 


J0rgensen & Hjorth (1997), 








J0rgensen et al. (1999) 


A665 


0.18 


6 


J0rgensen k. Hjorth (1997), 








J0rgensen et al. (1999) 


CL1358+62 


0.33 


10 


Kelson et al. (1997) 


MS1512+36 


0.37 


2 


Bender et al. (1998) 


A370 


0.37 


7 


Bender et al. (1998) 


CL0024-M6 


0.39 


8 


van Dokkum & Pranx (1996) 


MS2053-04 


0.58 


5 


Kelson et al. (1997) 


MS1054-03 


0.83 


8 


van Dokkum et al. (1998) 
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As examples of high redshift studies of the FP we show in Fig- 
ure 18 the FP for the Coma cluster and for five clusters with redshift 
larger than 0.1. The data for the Coma are from J0rgensen (1999) 
and J0rgensen et al. (1995). Abell 2218 and Abell 665 are discussed 
by J0rgensen et al. (1999). The sources for the rest of the clusters 
are given on the figure. 




Fig. 18. The FP edge-on for Coma, A2218, A665, CL1358+62, 
CL0024-hl6, and MS2053-04. The sources of the data are given on the 
panels ('This paper' refers to J0rgensen et al. 1999). The skeletal symbols 
on panel (c) and (d) are the E-l-A galaxies. The photometry is calibrated 
to Gunn r in the rest frames of the clusters. The mean surface brightness 
log(/)e = -0.4((;u)e - 26.4) is in units of Lo/pc^ (cf. Section 2.4.1). 
The photometry is not corrected for the dimming due to the expansion 
of the Universe. The effective radii are in kpc (Hq = 50kms~^ Mpc~^ 
and Qo = 0.5). The solid lines are the FPs with coefficients adopted 
from J0rgensen et al. (1996), and with zero points derived from the data 
presented in the figure. Typical error bars are given on the panels; the 
thin and thick error bars show the random and systematic uncertainties, 
respectively. (Prom J0rgensen et al. 1999.) 
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From the change in the zero point of the FP as a function of 
redshift J0rgensen et al. (1999), in agreement with other studies, 
find that the M/ L ratios of the E and SO galaxies change very slowly 
with redshift. The sample of clusters shown on Figure 18 span about 
half of the current age of the Universe (for go — 0.5). Under the 
assumption that the galaxies evolve passively over this time interval, 
e.g. no merging and no formation of new E and SO galaxies, then 
it is possible to put limits of the redshift at which the majority of 
the stars were formed. This redshift is called the formation redshift. 
The study by J0rgensen et al. (1999) as well as other studies conclude 
that the formation redshift is larger than about 2.5 for qo = 0.5 and 
larger than about 1.5 for qo = 0.15. 

It is important to keep in mind that the assumption regard- 
ing passive evolution represents a very simplified view of the galaxy 
evolution. Most likely the real evolution since z ^ 0.6 cannot be 
modeled with passive evolution. If the galaxies experience on-going 
star formation, then the observed evolution will appear smaller than 
for passive evolution because the galaxies are continuously forming 
young bright stars. A similar effect can be caused by a series of 
smaller bursts of star formation. Finally, the interactions, the pos- 
sible merging and the morphological evolution of the galaxies over 
the last half of the age of the Universe cannot be ignored. Some of 
the E and SO galaxies that we observe in low redshift clusters may 
not have ended up in the samples if we could have observed those 
clusters at a much earlier stage in their evolution, simply because 
some of the E and SO galaxies may have formed recently by merging 
of spiral galaxies. 

7. SUGGESTED FUTURE PROJECTS 

We end this paper by a brief summary of some of the projects 
that may be carried out building on the techniques and results dis- 
cussed in this paper. We concentrate on projects that involve pho- 
tometry only. Some of the projects may be carried out using existing 
archive data from HST. 

Evolution of morphology as a function of redshift 

Dressier et al. (1997) have recently used HST/WFPC2 data of 
clusters to study the morphology-density relation as a function of 
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redshift. Dressier et al. found that the fraction of SO galaxies is 
lower at high redshift than at low redshift, while the fraction of 
spirals is higher at high redshift than at low redshift. However, this 
study is based on the traditional method of classifying galaxies. We 
suggest that a quantitative approach is taken to any study of the 
morphology. Of special interest would be to study how the relative 
disk luminosities LoZ-^tot for E and SO galaxies evolve with redshift, 
both in clusters and in the field. The first study of this kind was done 
for the cluster CL0024+16 {z = 0.39) by Bergmann & J0rgensen 
(1999) who found that the Ln/L^oi distribution for the E and SO 
galaxies in CL0024+16 shows a paucity of disk-dominated galaxies 
when compared to the Coma cluster. 

It would be valuable to apply the same technique and derive 
L-o/Ltot for a larger sample of cluster and field galaxies at redshifts 
larger than 0.1 to establish the possible evolution of the distribution 
of Ld/L^o^. This may be done using the HST/WFPC2 archive data. 

Studies of global colors 

There have been many studies of the global colors of galaxies as 
a function of redshift (e.g. Stanford, Eisenhardt & Dickinson 1998; 
Bower, Kodama & Terlevich 1998; Kodama et al. 1998 and references 
in these papers). However, most studies concentrate on the optical 
colors of the galaxies, while the near-IR (JHK) data is very sparse. 
The study by Stanford et al. includes near-IR data and addresses the 
question of how the color-magnitude relations for the near-IR colors 
evolve with redshift. 

The combination of optical and optical-infrared colors may be 
used to break the age-metallicity degeneracy (see Section 5). Ob- 
servations of low redshift E and SO galaxies in the near-IR may be 
used to establish the zero redshift properties and the methods and 
models needed to break the age-metallicity degeneracy. For high 
redshift galaxies {z > 0.5), the near-IR photometry may be obtained 
with 8-meter class telescopes with superior spatial resolution. Such 
data will give the possibility of studying the mean ages and mean 
metallicities as functions of redshift for significantly fainter galaxies 
than it is currently possible by obtaining spectroscopy with 8-meter 
class telescopes. 
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Color gradients 

While we have not discussed color gradients in this paper, color 
gradients provide an alternative method of studying galaxy evolu- 
tion. The color gradients in E and SO galaxies reflect underlying 
radial gradients in the metallicity (and maybe the age) of the stel- 
lar populations. Models for galaxy formation predict the sizes of 
these gradients. In general, the predicted gradients are steeper for 
models based on a monolithic collapse (Carlberg 1984) than for mod- 
els based on the merger hypothesis (White 1980). Determination of 
color gradients for high redshift galaxies requires high signal-to-noise 
data with very good spatial resolution. Several of the rich galaxy 
clusters observed with HST/WFPC2 have sufficiently high signal- 
to-noise data that a study may be carried out using the available 
archive data. 
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